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ABSTRACT 

The Southern Coalsack is located in the interior of the Upper Centaurus-Lupus (UCL) super bubble 
and shows many traits that point to a much more energetic environment than might be expected from 
a dark, starless molecular cloud. A hot, X-ray emitting, envelope surrounds the cloud, it has a very 
strong internal magnetic field and its darkest core seems to be on astronomical time scales "just about" 
to start forming stars. In order to probe the magnetic environment of the cloud and to compare with 
the optical/near infrared polarimetry-based field estimates for the cloud, we have acquired Faraday 
Rotation measurements towards the pulsar PSR J1210— 6550, probing the magnetic field in the vicinity 
of the cloud, and a comparison target, PSR J1435— 5954, at a similar line of sight distance but several 
degrees from the cloud. Both lines of sight hence primarily probe the UCL super bubble. The earlier 
estimates of the magnetic field inside the Coalsack, using the Chandrasekhar-Fermi method on optical 
and near-infrared polarimetry, yield B^ = 64-93 jaG. However, even though PSR J1210— 6550 is located 
only —30 arc minutes from the (CO) edge of the cloud, the measured field strength is only By = 1.1±0.2 
fiG. While thus yielding a very high field contrast to the cloud we argue that this might be understood 
as due to the effects on the cloud by the super bubble. 

Subject headings: ISM: magnetic fields — ISM: individual: Southern Coalsack — pulsars: general — 
pulsars: individual (PSR J1047-6709, PSR J1210-6550, PSR J1435-5954) 



1. INTRODUCTION 

The Southern Coalsack provides one of the most spec- 
tacular sights on the southern night sky - at least for any 
interstellar medium astronomer - with its very dark neb- 
ulosity contrasted by the Southern Cross and the Jewell 
Box cluster (NGC 4755). Beyond its naked eye attrac- 
tiveness, it is a well-studied dark cloud, which, due to its 
relative vicinity (high spatial resolution) and location in 
the Galactic plane (abundance of background stars suit- 
able for optical and UV spectroscopy), provides one of 
the best test cases for the study of star-less molecular 
clouds (fo r a thorough introduction to the cloud see the 
review by Nvman (2008)). The cloud has been mapped in 
12 CO by Nvman_etaL| (11989|), who foun d a total mass of 
-3500 MO, and in 13 CO bv feato et aD ([I999h . The dis- 
tance to the cloud (d) has been discussed by a number o f 
authors, including ISeidensticker fc Schmidt- Kalerl (| 19891 ) 
who performed an extensive spectroscopic and photomet- 
ric survey of field stars in the region and find that the cloud 
is made up of two components, one at d — 188 pc and the 
other at d —243 pc. Later authors, generally, do not sup- 
port the dual components, but do find a cloud distance 
i n bro a d agreement with th e se authors (see e.g. iFrancol 
([1989ft : IStraizvs et al.l (Il994h: iKnude fc Hod (Il99^ ). The 
study by ISeidensticker fc Schmidt-Kalerl ([1989) had the 
additional advantage of mapping out the extinction well 
beyond the cloud and they conclude that the space beyond 
the Coalsack is relatively free of additional extinction until 
the Carina Arm is reached at —1.3 kpc. The Coalsack is 
located at /— 300°. W ith a nominal spiral arm pitch angle 
if -11° (Vallee 2005) the relative angle between the line 
of sight and the undisturbed magnetic field direction is of 
the order of 20-30°. 



ICrawfordl (Il99lh pointed o ut that, based on the H I ob- 
servations by Ide Geusl (|1992[ ). the cloud is located inside 
the Up per Centaurus-Lupus (UCL) super bubble. Ide Geusl 
(|l992f ) estimates that the super bubble has a radius of 
R=110±10 pc, an estimated age about 11 Myr, an es- 
timated total energy input of — 0.9xl0 51 ergs, and de- 
pending on the adopted distance to the Coalsack, the shell 
w ould have overtaken the cloud about 2-5 Myr ago. 

I Andersson et al.l (|2004l ) found that the extended enve- 
lope of the cloud contained the high-ionization state O VI 
i on, characteristic of a gas temperature of —300,000 K 
([Shapiro fc Mo or el 1 19761 ). Because O V has an ionization 
potential of —113 eV, O VI is generally expected to be 
produced by collisional ionization in the ISM. They also 
found that the cloud envelope could be seen in soft X-rays 
from the ROSAT PSPC ob s ervati ons. Based on these ob- 
servations I Andersson et al.l (|2004l ) proposed that the hot 
cloud envelope was d ue to the cloud ' s enve lopment by the 
UCL super bubble. [D uncan et aD (|l995f ) noted the ex- 
istence of a 2.4 GHz continuum emission arc around the 
cloud, which they d esignated G303.5+0, and the data from 
iDuncan et al.l (| 19971) show that this em ission is polarized 
(see Figure 2). IWalker fc Zealevl (|l998f ) identified this ra- 
dio feature with an Ha shell that they named "The Coal- 
s ack Loop" . 

lAndersson fc Potterl (|2005f ) used opti cal polarimetry 
and a Chandrasekhar-Fermi (CF) analysis ([Chandrasekhar , 
Il953l ) to estimate the plane-of-the-sky magnetic field 
strength in the cloud and found a surprisingly large field 
of B^=93±23 /iG, but one which is consistent, in equipar- 
tition, with the thermal p r essur e of the X-ray e mittin g gas 
seen by lAndersson et al.l ([2004D. lLada et al.l (|2004t ) used 
the H-band polarimetry of I Jones et al.l ([1984 ) and results 
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Fig. 1. — A sketch of the relevant Galactic feat ures is shown, w ith the Sun in the upper right. The UCL and the Coalsack are located 
based on the observat ions in lde Geua (Q.992) and Crawford (1991). The pulsar distances are based on the NE2001 electron d ensity mode l 
(Cordes & Lazio 2002). The gray area in the lower part of the figure represents the Car ina spiral arm with a pitch angl e of 11 ° ([Vallee 2005) 
and a dista nce, in the direction of the Coalsack, of 1.3 kpc based on the results from Seidensticker Sz Schmidt-Kaler (1989). The HI shell 
detected by McClure-Griffiths et al. (2001) is indicated by a dashed circle. 



from their own C 18 (J=2-l) observations to perform a CF 
analysis to estimate the magnetic field strength in Tapia's 
Globule 2, the darkest of the cores in the cloud. They 
estimate a polarization angle dispersion of <jq=0.66 rad 
(~38°) and hence derive a plane-of-the-sky magnetic field 
strength of ~24 jiG. We will argue below, in sec- 
tion 13.21 that separating the polarization data probing 
Tapia's Globule 2 from the surrounding cloud material, 
a narrower dispersion of polarization angles is found and 
hence a higher magnetic field estimate for the globule 
of B± ^65 /iG, consistent with the optical polarimetry 
es timate. 

lLada et ail (|2004f ) also used near-infrared photometry 
to map the density structure in Tapia's Globule 2. They 
detect a ring-like column density enhancement that they 
interp ret as a contrac t ing s hell of enhanced space den- 
sit y. IHennebelle et al.l (|2QQ6l ) modeled the structure seen 
by lLada et al.l ( 2004 ) and concluded that it is indeed due 
to a contracting shell, which they predict will lead to an 
onset of star formation about 10 4 years from now. They 
further show that the structure requires an external pres- 
sure transient to have passed over the cloud about 1.2 x 10 6 
years ago to trigger the cloud collapse. They speculate 
that the triggering event might have been the collision 
of two cloud fragments, giving rise to Tapia's Globule 
2. Models of cloud compression an d collapse for magne- 
tized clouds ([Li fc Nakamural l2002f ) can provide a simi- 
lar column density pattern, but would then require the 
magnetic field to be o riented close to the line of sight. 
iRathborne et al.l (l2009t), while n ot able to rule out the 
model of IHennebelle et ail ([20061 ), prefer a model where 
the structure originates as a merger of two turbulent flows. 

The location of the Coalsack relative to the UCL super 
bubble, the hot envelope, the strong magnetic field in the 
cloud and the presence of a possibly contracting core in 
Tapia's Globule 2, therefore seems to lead to a coherent 



picture of the Southern Coalsack as a cloud overtaken, 
compressed and heated by the envelopment into a super 
bubble and thence as a nearby example of star formation 
triggered over large distances. 

To further probe this scenario we have acquired Fara- 
day rotation measurements towards the pulsars PSR 
J1210-6550 and PSR JJ1435-5954. The pulsar PSR 
J1210— 6550 is located only about 30' from the lowest 
12 CO c ontour in the South- W esterly- most cloudlet in the 
map of iNvman et~ai1 §989). The pulsar is at an esti- 
mated distance of D= l.lb kpc, based on t he NE2001 
electron density model ([Cordes fc Lazioll2002f ) and hence 
on the near-side of the Carina arm and should primar- 
ily ([Seidensticker fc Schmidt- Kalerl 1 19891 ) be probing the 
material associated with the UCL super bubble and the 
Coalsack envelope. For comparison, we also observed 
the pulsar PSR J1435— 5954, located about 17° further 
along the Galactic plane, but at a similar distance of 
D=1.1S kpc. Figure [1] shows the sketch of the geometry of 
the objects. Of course, it is important to remember the in- 
herent dichotomy between the magnetic field components 
probed by polarimetry and the CF analysis on the one 
hand (plane-of-the-sky component) and Faraday rotation 
measurements (line-of-sight component) on the other. 

Both these pulsars were am ongst later discove ries by 
the Parkes multibeam survey (Hob bs et al.ll2004D . PSR 
J1210— 6550 is a very faint, long-period (4.3 s) pulsar with 
a very small duty cycle (~0.1%), with a period- aver aged 
flux density at 1400 MHz (S1400 ) of ~ 0.1 mJy, whereas 
PSR J1435— 5954 is a relatively bright pulsar. Little is 
known about their properties apart from what has been 
deduced from their discovery and initial timing observa- 
tions. Additionally, we observed PSR J1047— 6709 in each 
observing session, a bright pulsar with a well-known rota- 
tion measure (RM), to serve as a 'control pulsar' for vali- 
dating polarimetric calibration and the RM determination. 
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Fig. 2. — Overlaid on a false-colour X-ray map of the Coalsack region (red=l/4 keV, green=3/4 keV, blue=1.5 keV; all ROSAT data) are: (i) 
the lowest detected CO (J=l-0) contour for the Coalsack (white line; Nyman, Bronfman Sz Thaddeus, 1989), and (ii) contours of the polarized 
2.4 GHz emission (Duncan et al., 1997; gray contours^ 50, 70, 90 and 110 mj per beam). The 2.4 GHz survey covers \b\ < 5° and we have 
here cut out data for 297° < I < 307°, producing a 10° x 10° area indicated by the dashed gray line. The locations of PSR J1210— 6550 and 
PSR J1435— 5954 are shown as star symbols. As shown by Andersson et al. (2004), the enhanced X-ray emission surrounding the Coalsack is 
due to a hot halo enveloping the cloud, which is most likely due to an interaction between the cloud and the Upper Centaurus-Lupus super 
bubble. The polarized radio emission points to a synchrotron origin, consistent with a hot magnetized plasma. 



The basic parameters of all three pulsars are listed in Table 
1, along with the measurements of RM and the estimates 
of magnetic field strength derived from our data. We note 
that even in the best available pulsar surveys, no further 
observable targets are available meeting our requirements 
of 1) on-the-sky proximity to the Coalsack and 2) a dis- 
tance less than that of the next spiral arm (Carina). 

The locations of PSRs J1210-6550 and J1435-5954 on 
the sky and their relation to the known features around 
the Coalsack region are further illustrated in Fig. [2j 
which shows an X-ray map in that area (from ROSAT 
data) along with t he detections of CO (J=l-0) from 
iNvman et al.l (Il989l ) and the polarized 2.4 GHz emission 
from lDuncan et al.l ([1997D . 

1.1. Magnetic Fields from Rotation Measure 
Observations 

Rotation measure quantifies the degree of Faraday rota- 
tion that electromagnetic waves undergo as they propagate 
through the ISM from the pulsar to the Earth. This ro- 
tation is caused by the interaction of the electromagnetic 
waves with the magnetized plasma of the ISM and more 
specifically along the line of sight to the pulsar. For a pul- 
sar located at a distance D, the rotation measure (RM) is 
given by 

RM=^f D n em (l).dl (1) 



where dl is the path vector element in the direction of wave 
propagation; e and m e are the charge and mass of electron, 
respectively; n e is the free electron density, B is the mag- 
netic field vector, and c is the speed of light. The above 
equation can be simplified to 



RM = 0.812 / n e (l)B(l).d\ 
Jo 



(2) 



where D is in parsecs, n e is in per cubic centimeter, B is in 
microgauss and RM is in radians per square meter. The 
integral of n e along the Earth-pulsar line of sight is the 
dispersion measure (DM), and is given by 



DM 



/ n e 
Jo 



(l)dl 



(3) 



and is usually expressed in units of parsecs per cubic cen- 
timeter (pc cm -3 ; i.e. D in pc and n e in cm -3 ). 

The quantity DM can easily be determined from the ar- 
rival time delays measured at different frequencies across 
the observing observing band (e.g. Lorimer & Kramer 
2005) and, consequently, for all catalogued pulsars which 
are observable in radio this is a well-known quantity. Thus, 
a measurement of RM enables a direct estimation of the 
magnetic field strength weighted by the free electron den- 
sity. The mean value of the line-of-sight component of B is 
thus given by 



(B||> 



JO n * 



(l)B(l).dl 



{l)dl 



1.232 



RM 
~DM 



(4) 
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Table 1 

Basic parameters and measured quantities of observed pulsars 



Pulsar 


PSR J1047-6709 


PSR J1210-6550 


PSR J1435-5954 


Pulse Period (ms) a 


198.4514 


4237.0102 


472.9954 


Dispersion Measure (pc cm -3 ) a 


116.156 ±0.002 


37 ±6 


44.26 ±0.11 


Distance (kpc) 6 


2.88 


1.15 


1.18 


Galactic longitude (°) 


291.31 


298.77 


315.58 


Galactic latitude (°) 


-7.13 


-3.29 


0.39 


Flux at 1400 MHz S1400 (mJy) 


2.2 ±0.3 


0.2 ±0.1 


1.15 ±0.35 


Linear polarization (%) 


75 


35 


10 


Pulse width W50 (nis) 


2.1 ±0.1 


35 ±1 


18 ± 1 


Rotation measure (rad m -2 ) 


-84 ±3 


-32 ±4 


-24 ±7 


Magnetic field strength B\\ (/iG) 


-0.95 ±0.03 


-1.1 ±0.2 


-0.7 ±0.2 



2 From the ATNF Pulsar Catalog (also s ee iHp^bseTaTl (12004 ) 
6 Estimate based on the NE2001 model ([Cordes fc Laz io 2002) 



The rest of this paper is organized as follows. We briefly 
describe our observations and data reduction in § [2J deter- 
mination of rotation measures in ? 13. H and a re-analysis of 
near infrared polarization data in § 13.21 In the later sec- 
tions (§H]and[5j) we present our main results and discuss 
their implications for the related physical models. Our 
conclusions are presented in § [6j 

2. OBSERVATIONS AND DATA REDUCTION 

Pulsar observations were made us ing the central beam 
of the 13-beam multibeam receiver ([Stavelev- Smith et al.l 
1996) on the Parkes 64- metre telescope. This receiver is a 
dual-channel cryogenic system sensitive to orthogonal lin- 
ear polarizations with a system equivalent flux density of 
36 Jy. In all observing sessions, signals centered at 1369 
MHz with a bandwidth of 256 MHz were recorded using 
the Parkes Digital FilterBank 3, which is capable of pro- 
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Fig. 3. Polarimetric profiles of PSR J1047-6709 at 1.4 GHz 
from observations made with the Parkes multibeam receiver (MJD 
= 55190). The blue and red lines represent circular and linear po- 
larization respectively, while the black lines represent the total in- 
tensity. Only 80° of pulse rotation is shown, and absolute position 
angles at infinite frequency are shown after correction for a rotation 
measure of —80 rad m -2 . 



ducing high-resolution profiles in both time and frequency 
(up to 1024 phase bins and 2048 spectral channels). For 
all our observations, the full Stokes profiles were recorded 
across the 256 MHz bandwidth split into 1024 frequency 
channels, with a 512-bin resolution across the pulse period, 
over 10-second sub- integrations. The combination of such 
short sub-integration and high time and frequency resolu- 
tions enable a high efficiency in excising the data segments 
corrupted by radio frequency interference (RFI). 

We collected data through multiple observing ses- 
sions conducted between December 2009 and April 2010. 
The integration times varied depending on the antici- 
pated flux densities of pulsars; PSR J1210— 6550 was typ- 
ically observed for durations of 1-2 hr, whereas PSR 
J1435-5954for 0.5-1 hr and PSR J1047-6709for 15-30 
min. Due to its low flux density PSR J1210— 6550 requires 
long integrations to enable high quality detections (~ 1 
hr to reach S/N ~ 40). Each pulsar observation scan was 
preceded by a short (2-min) observation of the pulsed cal- 
ibration signal, which is linearly polarised and broadband, 
and injected into the feed at 45° to the two signal probes. 
In addition, we made observations of the radio galaxy Hy- 
dra A (3C218), which is assumed to have a flux density of 
43.1 Jy at 1.4 GHz, in each observing session in order to 
enable flux calibration of the data. 

For off-line data processing, we u sed the psrchive pul- 
sar data analysis software package (Hot an et al.ll2004} ). In 
brief, the raw data were first subjected to a process of 
RFI excision, whereby any frequency channels containing 
strong interference and sub-integrations affected by strong 
impulsive interference were given zero weight. In addition, 
the upper and lower 5% of the band at the band edges were 
given zero weight due to a low gain and other instrumental 
effects in these ranges. Using the pulsed calibration sig- 
nal as reference, variations in instrumental gain and phase 
across the band were removed and corrected for parallac- 
tic angle variations. The data were then summed in time 
to obtain four integrated Stokes profiles for each frequency 
channel. A full polarimetri c calibration m odeling of the re- 
ceiver as discussed in Ivan Stratenl (|2004l ) is not critical for 
the RM determination and therefore not considered in our 
analysis. The position angles were then computed for each 
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phase bin from the Stokes Q and U parameters and the 
flux-density scale was determined from the observations of 
the Hydra calibration observations. 

3. ANALYSIS 

3.1. Determination of Rotation Measures 

Polarimetric profiles obtained from our analysis are 
shown in Figs. [3] and H While PSR J1210-6550is very 
faint, it shows a substantial degree of linear polarization 
(-35% of the total intensity). PSR J1435-5954is rela- 
tively bright in comparison, however the degree of polar- 
ization is much lower (linear and circular polarizations are 
—8% and —3% of the total intensity, respectively). The 
signal-to- noise ratios (S/N) of the linearly polarized in- 
tensities are therefore not very large and limit the signif- 
icance achievable in the RM estimates. Flux densities of 
both pulsars show significant variations with time, pre- 
sumably due to long-term interstellar scintillation effect s 
expected at such moderate DMs (e.g. Bh at et all ll999), 
and the quoted numbers are mean values from multiple 
measurements made over time spans of several months. 

Our initial estimates of RM were determined using the 
standard procedure, where we searched for a peak in the 
total linearly polarized intensity L = (Q 2 + U 2 ) 1 ^ 2 ob- 
tained by summ ing the calibrated data in frequency (e.g. 
lHan et a l. 2006). This search is carried out over a large 
range of RM, up to ± 1000 rad m -2 , in steps of 1 rad m -2 . 
Further, using the RM value corresponding to the peak, we 
summed the data to form the upper and lower band pro- 
files. These were then used to refine the RM estimate by 
taking the weighted mean of the position-angle differences 
between the upper and lower bands across the profile, with 
the weight inversely proportional to the square of the error 
in position-angle difference for each pulse phase bin. The 
RM is then recomputed and the procedure is iterated until 
convergence. 

In order to confirm and cross-check our initial esti- 
mates, as well as to obtain more robust RM measurements, 
we applied th e RM determination method developed by 
iNoutsos et al.l (|2008l ), where a quadratic function was fit- 
ted to the PA vs frequency across the full observation 
band. These position angles are computed from the Stokes 
parameters Q and U that are summed across the phase 
bins corresponding to the "on" pulse. This quadratic fit- 
ting algorithm finds the best fit by means of a Bayesian 
likelihood test, and uses a fitting function of the form, 
PA = PA + c 2 RM(l/f 2 - 1//2), where PA is the PA 
at frequency fo and fj is the frequency channel j. The 
main advantage of this method is that it accounts for the 
180° ambiguity inherent in the computation of the mean 
PAs from the Stokes Q and U. The free parameters PAq 
and RM were stepped through the ranges (0, n) rad and 
(-1000, 1000) rad m ~ 2 in 1° and 1 r ad m~ 2 , respectively. 
As demonstrated bv lNoutsos et al.l (|2008f h this method is 
more robust and yields more reliable RM estimates. The 
main caveat however is that the signal-to- noise ratio (S/N) 
per frequency channel or sub-band needs to be sufficiently 
large for this method to yield meaningful results. While 
we were able to apply this method successfully to both 
PSRs J1047-6709 and J1210-6550 (cf. Fig 5), it was not 
feasible in the case of PSR J1435— 5954 owing to the very 
low signal-to- noise of its linearly polarised flux (cf. Fig 4). 
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Fig. 4. Polarimetric profiles of PSR J1210-6550 (top) and PSR 
J1435— 5954 (bottom) at 1.4 GHz from observations made with the 
Parkes multibeam receiver (MJD = 55190 and 55302 respectively). 
The blue and red lines represent circular and linear polarization 
respectively, while the black lines represent the total intensity. PSR 
J1210-6550 is moderately polarized, whereas for PSR J1435-5954 , 
the measured linear polarization is only ~8% of the total intensity. 

As a result, the RM estimate of this pulsar is derived from 
the standard method as discussed earlier in this section. 

The final estimates of RM derived from our analy- 
sis are summarized in Table 1. For the control pulsar 
PSR J1047— 6709, we obtained four independent mea- 
surements over a time span of four months, which are 
found to be consistent within their 1 a measurement 
uncertainties, with a mean value of — 84 ± 3 rad m -2 . 
For comparison, values from the published literature 
are -73 ± 3 rad m~ 2 dHan et al.l l2QQ6h and -79 ± 2 
rad m~ 2 (|Noutsos et al.l 120081 ). Thus, our RM estimate 
is consi stent (within 2 a) w ith the more recent measure- 
ment of INoutsos et al.l ([2008), though somewhat larger. 

For both PSR J1210-6550and PSR J1435-5954, we 
have obtained 2-3 independent measurements (with high 
significance) over a time span of 3-4 months. The re- 
sultant average values are — 32 ± 4 rad m -2 and — 24 ± 7 
rad m- 2 for PSR J1210-6550and PSR J1435-5954, re- 
spectively. The relatively low significance on the RM 
estimate of the latter is primarily to do with its low 
degree of polarization (~8%). In comparison, PSR 
J1210— 6550 exhibits a moderately high linear polariza- 
tion (~35%) and hence its RM estimated with a relatively 
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higher significance (Fig. [5]). However, due to its very nar- 
row pulse (duty cycle of ~1%), typically only four phase 
bins span the on pulse emission with our 512-bin reso- 
lution across the pulse period. The consistency between 
multiple measurements of each pulsar and the general 
agreement seen between the published and our measured 
RM estimates for the control pulsar give us confidence on 
the reliability of our RM measurements. 

For the current analysis, the dispersion measures to- 
wards each of the pulsars under consideration were take n 
from the ATNF Pulsar Catalog ([Manchester et al.ll2005l ): 
for both PSR J1210-6550and PSR J1435-5954, the DM 
measurements (37 ±6 p c cm" 3 and 44.26 ± 0.011 pc cm -3 , 
respectively) are from iHobbs et al.l (|2004l ) that reported 
the original discoveries of these objects. 

3.2. The Magnetic Field in Tapia's Globule 2 
As discussed in the introductio n, lLada et al.l ( 2004) used 



near infrared polarimetry from I Jones et al.l (|1984f ) and 
their own deep 18 CO observations of Tapia's Globule 2 
to estimate the magnetic field strength in this dark core. 
We here discuss this r esult in the light of a reanalysis of 
the polarization data. I Jones et"aD (|l984h note about their 
H-band polarization map that: 

"Immediately apparent in Figure 4 is a definite change 
in polarization angle between the Northern block (block E 
[...]) and the block centered on the globule (block D). Also 
the polarization within the 120' radius circle [the globules 
"half- extinction radius"] is about a factor 2 higher than in 
block E". 

Fig ure |6] shows the distribution of polarization angles 
from I Jones et al.l (|l984f ) where we have color coded the 
data from "block D" in blue (dark grey) and those from 
blocks E and F in yellow (light grey). The change in the av- 
erage value of the polarization angles between the areas is 
clearly seen. A student's t-test analysis comparing the po- 
larization angles of the "block D" sample vs. the combined 
"block E" and "block F" samples yields a t-probability of 
P=0.0003 and hence shows that the two populations are 
distinct. The position angle dispersion for the full s ample 
is cr#=39±10°. However, if we follow Uones et al.l (|1984[ ) 
and separate out only the "block D" polarization as trac- 
ing the field in Tapia's Globule 2, we derive a position 
angle dispersion of gq =14±2°. 

4. RESULTS 

The estimates of By derived from our RM measure- 
ments are summarized in Table 1. The measured field 
strength along the sightline toward PSR J1210— 6550 is 
only — 1.1 ± 0.2 /iG , which is much lower than the plane 
of the sky component (B_l = 93 ± 23 /iG ) estimated 



by lAndersson fc Potterl ((2005). The field strength esti- 
mated toward PSR J1435— 5954 is somewhat lower, By = 
— 0.7 ± 0.2 /iG , but again is directed away from us. 

For our reanalysis of the NIR polarimetry in Tapia's 
Globule 2, if we retain the values of turbulent velocity 
(a v =0.23 km s" 1 ) an d space density (n=10 4 cm -3 ) derived 
bv lLada et all (|2004f h but use the narrower position angle 
dispersion, we estimate a plane-of-the-sky magnetic field 
strength for Tapia's Globule 2 of B^ ~ 64 /iG. This higher 
value is consistent (within 1. 5cr) with the field strength 
(B ± =93±23 /iG) derived bv lAndersson fc Potterl (2005) 
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Fig. 5. — Position angle vs frequency plots for PSR 
J1047-6709and PSR J1210-6550 , illustrating the RM estimation 
using the quadratic fit method (see § 13.11 for details). The best 
fit RM determined toward PSR J1210-6550is -32 ± 4 rad m~ 2 , 
which translates to By ~ — 1.1 zb 0.2 /iG . These position angles are 
computed from the Stokes Q and U. For PSR J1210-6550 , they 
are summed over the four on pulse phase bins that show polarized 
intensity above a 5-cr threshold. The large errors on the PAs of this 
pulsar are thus due to its relatively low degree of polarization and 
a small number of on pulse phase bins. 



using optical polarimetry of the full Coalsack cloud. It is, 
however, si gnificantly higher than a value of & 25 /iG 
derived by lLada et all (|2QQ4f ) based on equipartition ar- 
guments between magnetic, gravitational and kinetic en- 
ergies, which is also very close to their CF-analysis based 

v alue. 

iRathborne et al.l (|2QQ9h have used a detailed analysis 
of the C 18 (J=l-0) emission to argue that the globule 
is in fact the confluence of two subsonic flows. They 
estimate the average space density of the globule to 
n ^2.7xl0 3 cm -3 . They, however, also detect CS emis- 
sion from the inner parts of the globule. The CS molecule 
has a high critical density and is expected to trace high- 
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density gas (n >10 5 cm -3 ). Finally, the spatial resolution 
of Jones et al's polarimetry data is not high enough to 
fully sample the m agnetic field variations on the scale of 
the observations of iRathborne et all (|2009f ). Hence, a de- 
tailed understanding of the magnetic field in the core of 
the globule will require additional observations. 

5. DISCUSSION 

A number of theoretical studies have addressed the char- 
acteristics of super bubble evolu tion ([Weaver et al.|[l977t 
iTomisakal 119981 : IStil et al.l 120091 ) and of the effects of a 
magnetized interstella r cloud being overtaken by a super- 
nova remnant (SNR) (jMac Low et al.l 1 19941 : iMelioli et all 
2005; Leao et al. 2009). While there are no detailed mod- 
eling studies specifically addressing a magnetized cloud en- 
veloped in the relatively slowly expanding super bubble, 
we can gain some insights into the physics, and estimate 
the characteristics, of such a cloud by viewing it as a de- 
tached part of the inner surface of a super bubble wall 
or by comparing to the studies of a cloud enveloped by a 
SNR. As shown by m ost super bubble models starting with 
I Weaver et al.l (|1977[ ), the inside of the super bubble shell 
will reach temperatures in excess of 10 6 K and should reac h 
thermal pressures of P^IO 6 K cm -3 ([Tomisakal 1 19981 ). 
While the magnetic field is swept up and enhanced, only 
marginal magnification of the field strength i s exp ected 
for the in side of the bu b ble s hell ([Tomisakal Il998f ) . In 
contrast, iMac Low et al.l ([19941 ) find that the field comes 
into equipartition with the gas pressure for the cloud over- 
taken by a SNR and that the magnetic field prote c ts the 
cloud from being shredded by the shock. iLeao et al.1 (2009) 
have calculated the star formation efficiency for clouds en- 
veloped by a SNR driven by an energy injection of 10 51 erg, 
but only briefly discuss the results for SNR in the radia- 
tive phase, where the expansion velocity is more compara- 
ble to that for the UCL super bubble. Nonetheless, based 
on their Figure 9, the pre dicted imminent onset of star 
formation in the Coalsack (He nnebelle et al1l2006f ) is con- 
sistent with being triggered by the interaction with the 
UCL super bu bble. 

As shown bv ITomisakal (|l998l ) and lStil et all (|2009f ). an 
expanding super bubble will sweep up the ambient mag- 
netic field, compressing it in the direction perpendicular to 
the undisturbed field direction and leaving an almost field- 
free bubble interior. Those calculations show a precipitous 
drop of the magnetic field at the inside of the bubble wall. 
Hence, by treating the cloud as a detached part of the wall, 
we expect that the magnetic field frozen in to the cloud will 
also sh ow an abrupt drop-off. As noted bv lMac Low et all 
(1994) the initial compression of a magnetized cloud is ex- 
pected to reverse on time scales short compared to the 
lifetime of the bubble. Moreover, in the simulations of 
super bubbles, the magnetic field leads to relatively thick 
shells - partic ularly in the di rection per pendicular to the 
ambient field (IStil et al.ll2009h . However. IStil et al.l (|2009f ) 
point out that including realistic cooling functions leads 
to thinner shells and stronger magnetic fields than in adi- 
abatic simulations. Since the average magnetic field, as 
traced by the optical polarimetry, runs NE-SW, such cor- 
rections would likely be particularly important for the lo- 
cation of PSR J1210— 6550, where cooling flows from the 
dense cloud likely produces a much more efficient cooling 
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Fig. 6. — The distribution of H-band polarization angles from 
Uones et al.l (|1984l ) shows a distinctly bi-modal distribution between 
the stars behind Tapia's Glubule 2 ( "block D" ) and the surrounding 
Coalsack material. For the full sample a position angle dispersion 
of A6 ~39° is deduced. If only the stars in "block D" are used, an 
angle dispersion of AO ~14° is found implying a magnetic field of 
B ± ^65 iiG. 



than in the locations where the cloud and surrounding hot 
gas are connected across the magnetic field lines. 

Because the field inside the super bubble in the 
IStil et all (|2QQ9f ) simulations is weak and turbulent, most 
of the predicted r otation measure will be due to the shell. 
On a large scale, IStil et all (|2QQ9l ) therefore estimate the 
rotation measure for a 10 Myr old bubble, seen for sight- 
lines perpendicular to the undisturbed field, to be less 
than ~ 20 rad m -2 , dominated by asymmetries in the 
shell. For sightlines along the undisturbed field, relatively 
large rotation measures are seen close to the tangent of 
the shell while the rotation measure in the projected cen- 
ter of the bubble significantly lower. These predictions 
are generally consisten t with our observation s through 
the UCL super bubble. IVallee fc Bignelll (|l983f ) observed 
extragalactic background sources through the Gum Neb- 
ula and derive an average RM in the inner part of the 
Gum Nebula of ~130 rad m -2 . However, the Gum Neb- 
ula is located at / ~270° and hence much closer to the 
pr ojected direction of the Galactic magnetic field. Also, 
as IVallee fc Bignelll (pj)83) used extragalactic background 
sources, contributions from the rest of the Galaxy along 
th e line of sight ca nnot be excluded. As the simulations 
bv lStil et al.l (|2009f ) show, for this geometry a fairly large 
RM is expected. 

The line of sight towards PSR J1435-5954 is located 
at a larger Galactic longitude (7^316°) than that of PSR 
J1210— 6550 and hence at a larger angle relative to the 
nominal undisturbed Galactic magnetic field. The lower 
RM seen towards this pulsar, as compared to the one for 
PSR J1210-6550, is thus also in line with the models by 
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iStil et all (120091) . 

Based on the results from the CF a nalyses from optica l 
(1 Andersson fcPotter|[2005l) and NIR ([ Jones et all ([1984D : 
lLada et all (|2004f ); our Section 1X2]) , the results presented 
here thus seem to indicate that both the sightlines probed 
fall in the "diffuse" part of the super bubble, including the 
one towards PSR J1210— 6550. This then suggests that 
the edge of the magnetized Coalsack cloud is quite sharp 
and located close to the lowest CO contour - at least in 
the direction of the magnetic field. This is consistent with 
models, however needs to be tested in detail. Further ob- 
servations and models specifically targeted at simulating 
a magnetized cloud enveloped in a super bubble with pa- 
rameters tailored to those of the Coalsack and the UCL 
super bubble are required to address the viability of this 
interpretation. 

6. CONCLUSIONS 

We have acquired rotation measure data for two pulsars 
behind the Upper Centaurus-Lupus super bubble, one of 
which probes a sightline very close to the Southern Coal- 
sack. While constituting a small sample, these two are the 
only readily observable pulsars to provide a specific probe 
of the Faraday rotation in the 3-dimensional vicinity of 
the cloud. We find that the results for both lines of sight 
are consistent with the predictions of models for a magne- 
tized super bubble without internal clouds. Since earlier 
estimates of the magnetic field strength in the Coalsack in- 
dicate a strong field (at least in the plane of the sky), this 
indicates that the magnetized cloud is either very sharply 
bounded or that the field in the cloud is oriented almost 
completely in the plane of the sky. 

We have shown that the observational data set avail- 
able for the Coalsack is consistent with the hypothesis of 
a cloud relatively recently overtaken by a super bubble and 
therefore at the brink of star formation. Specific model- 
ing of the interaction of the Coalsack and the UCL will 
be required to confirm this picture; the proximity of the 
cloud and the wealth of observational data are promising 
to provide important constraints on such models. 
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